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fi]  The  variation  of  seaward  and  shoreward  infragravity  energy  fluxes  across  the  shoaling 
and  surf  zones  of  a  gently  sloping  sandy  beach  is  estimated  from  field  observations  and 
related  to  forcing  by  groups  of  sea  and  swell,  dissipation,  and  shoreline  reflection.  Data 
from  collocated  pressure  and  velocity  sensors  deployed  between  1  and  6  m  water  depth  are 
combined,  using  the  assumption  of  cross-shore  propagation,  to  decompose  the  infragravity 
wave  field  into  shoreward  and  seaward  propagating  components.  Seaward  of  the  surf 
zone,  shoreward  propagating  infragravity  waves  are  amplified  by  nonlinear  interactions 
with  groups  of  sea  and  swell,  and  the  shoreward  infragravity  energy  flux  increases  in  the 
onshore  direction.  In  the  surf  zone,  nonlinear  phase  coupling  between  infragravity  waves 
and  groups  of  sea  and  swell  decreases,  as  does  the  shoreward  infragravity  energy  flux, 
consistent  with  the  cessation  of  nonlinear  forcing  and  the  increased  importance  of 
infragravity  wave  dissipation.  Seaward  propagating  infragravity  waves  are  not  phase 
coupled  to  incident  wave  groups,  and  their  energy  levels  suggest  strong  infragravity  wave 
reflection  near  the  shoreline.  The  cross-shore  variation  of  the  seaward  energy  flux  is 
weaker  than  that  of  the  shoreward  flux,  resulting  in  cross-shore  variation  of  the  squared 
infragravity  reflection  coefficient  (ratio  of  seaward  to  shoreward  energy  flux)  between 
about  0.4  and  1.5.  INDEX  TERMS:  4546  Oceanography:  Physical:  Nearshore  processes;  4560 
Oceanography:  Physical:  Surface  waves  and  tides  (1255);  4594  Oceanography:  Physical:  Instruments  and 
techniques;  KEYWORDS:  infragravity  waves,  energy  flux,  surf  zone,  nonlinear  phase  coupling,  shoreline 
reflection,  cross-shore  evolution 


1.  Introduction 

[2]  Over  a  distance  of  just  a  few  wavelengths  on  a 
moderately  sloped  beach,  shoreward  propagating  sea  and 
swell  (periods  of  roughly  4-20  s)  steepen,  pitch  forward, 
break,  and  form  dissipative  bores.  This  evolution  is  accom¬ 
panied  by  the  generation  of  gravity  waves  with  periods  of  a 
few  minutes  (inffagravity  waves).  Here  field  observations 
are  used  to  estimate  the  variation  of  seaward  and  shoreward 
infragravity  energy  fluxes  across  the  shoaling  and  surf  zones. 

[3]  Previous  field  observations  have  shown  that  standing 
waves  can  be  a  dominant  component  of  the  nearshore 
inffagravity  wave  field  [e.g.,  Suhayda,  1974;  Huntley, 
1976;  Huntley  et  al.,  1981;  Oltman-Shay  and  Howd, 
1993].  Standing  inffagravity  waves  are  consistent  with 
relatively  weak  forcing  and  damping  and  strong  shoreline 
reflection.  However,  there  also  is  evidence  that  a  substantial 
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progressive  infragravity  component  can  result  from  non¬ 
linear  forcing  of  shoreward  propagating  infragravity  waves 
by  groups  of  sea  and  swell  [e.g,,  Munk,  1949;  Longuet- 
Higgins  and  Stewart ,  1962;  Meadows  et  al ,  1982;  Elgar 
and  Guza ,  1985;  List ,  1992;  Herbers  et  al,  1994].  Non- 
linearly  forced  shoreward  propagating  and  free  seaward 
propagating  (resulting  from  shoreline  reflection)  infragrav¬ 
ity  waves  are  expected  to  have  different  cross-shore  ampli¬ 
tude  variations,  leading  to  partially  standing  waves. 
Observations  in  13  m  water  depth,  2  km  offshore  of  Duck, 
North  Carolina  (the  site  of  the  observations  presented  here), 
show  that  the  net  (integrated  over  the  infragravity  frequency 
band)  energy  flux  can  be  directed  either  seaward  or  shore¬ 
ward  and  that  the  ratio  of  seaward  to  shoreward  fluxes 
(hereinafter  the  bulk  infragravity  reflection  coefficient  R2) 
deviates  substantially  from  1,  with  values  usually  between 
0.5  and  4.0  [Elgar  et  al,  1994;  Herbers  et  al,  1995b]. 
Reflection  coefficients  <1  (>1)  imply  an  energy  sink 
(source)  between  the  shoreline  and  13  m  depth  observation 
sites.  Three  R 2  regimes  were  identified.  With  low  sea-swell 
energy  the  infragravity  wave  field  was  dominated  by 
arrivals  from  remote  sources  that  were  either  partially 
dissipated  or  scattered  into  trapped  waves  shoreward  of 
13  m  depth,  and  R2  <  1.  With  moderate  sea-swell  energy, 
infragravity  energy  generated  onshore  of  13  m  depth  and 
radiated  seaward  exceeded  dissipation  and  trapping  losses, 
and  R2  >  1.  With  the  most  energetic  sea  and  swell  waves  the 
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13  m  depth  observations  were  near  the  surf  zone  where 
dissipation  might  be  important,  and  R2  <  1.  Nonlinear  phase 
coupling  of  infragravity  waves  to  groups  of  sea  and  swell 
seaward  of  the  surf  zone  (depths  between  about  8  and  200  m) 
has  been  investigated  with  observations  in  the  Atlantic  and 
Pacific  Oceans  [Herbers  et  al .,  1994,  1995b].  The  observed 
levels  of  phase  coupling,  estimated  with  bispectral  analysis, 
agreed  well  with  the  theory  for  second-order  bound  waves 
driven  by  quadratic,  nonresonant,  sea- swell  interactions 
[ Hasselmann ,  1962],  and  both  the  bound  wave  energy 
levels  and  the  ratio  of  bound  to  total  (free  plus  bound) 
infragravity  energy  increased  with  increasing  sea-swell 
energy  and  decreasing  depth.  Similar  increases  in  nonlinear 
phase  coupling  were  observed  in  pressure  measurements 
seaward  of  the  surf  zone  (in  3-6  m  depth)  [Ruessink,  1998]. 
However,  the  coupling  decreased  to  approximately  zero 
when  the  offshore  sea-swell  energy  was  so  large  that  the 
observations  were  within  the  surf  zone. 

[4]  Here  a  more  detailed  analysis  of  the  transformation  of 
infragravity  waves  through  the  shoaling  and  surf  zones  is 
presented  using  extensive  new  observations.  Collocated 
current  meters  and  pressure  sensors  deployed  in  1-6  m 
depth  (described  in  section  2)  are  used  to  decompose  the 
infragravity  wave  field  into  shoreward  and  seaward  prop¬ 
agating  components.  Errors  introduced  by  neglecting  the 
directionality  of  the  infragravity  wave  field  are  shown  to  be 
small  in  Appendix  A.  The  spatial  variation  of  infragravity 
energy  fluxes  is  related  to  forcing  by  groups  of  sea  and 
swell,  dissipation,  and  shoreline  reflection  in  section  3.  The 
directionality  of  infragravity  waves  is  discussed  briefly  in 
section  4.  A  more  detailed  analysis  is  presented  by  A. 
Sheremet  et  al.  (manuscript  in  preparation,  2002).  Results 
are  summarized  in  section  5. 

2.  Field  Experiment  and  Analysis 

[5]  Observations  were  obtained  from  August  to  Decem¬ 
ber  1997  during  the  Sandyduck  experiment  conducted  on  a 
sandy  beach  near  Duck,  North  Carolina.  A  two-dimensional 
array  of  electromagnetic  current  meters,  downward  looking 
sonar  altimeters,  and  pressure  sensors  was  deployed  in  1-6 
m  depth,  50-350  m  from  the  shoreline  (Figure  1).  Bathy¬ 
metric  surveys  obtained  a  few  times  a  week  with  an 
amphibious  vehicle  were  supplemented  with  nearly  contin¬ 
uous  altimeter  observations.  There  was  a  50  cm  high 
sandbar  crest  in  about  3.5  m  mean  water  depth  (cross-shore 
location  x  =  320  m)  and  a  transient  sandbar  in  about  1  m 
depth  (x  —  160  m)  (Figure  1).  Alongshore  variability  of  the 
bathymetry  usually  was  weak  over  the  instrumented  area, 
except  near  the  shoreline.  Redder  sen  et  al.  [2000]  and  Elgar 
et  al.  [2001]  give  additional  details  and  describe  conditions 
during  the  experiment. 

[6]  To  investigate  cross-shore  energy  fluxes,  surface 
elevation  time  series  of  shoreward  and  seaward  propagating 
infragravity  waves  (t/  and  r\~,  respectively)  are  constructed 
from  collocated  pressure  p  and  cross-shore  velocity  u  time 
series  using  the  assumptions  of  shallow  water  and  cross¬ 
shore  propagation  [Guza  et  al ,  1984;  Elgar  and  Guza , 
1985;  List ,  1992], 

=  ^-{p±u-Jhfg)  , 


Figure  1.  (top)  Plan  view  of  instrument  array.  Circles  are 
collocated  pressure  ( p )  and  horizontal  velocity  («  and  v) 
sensors.  Current  meters  usually  were  positioned  between  50 
and  100  cm  above  the  seafloor.  The  Field  Research  Facility 
coordinate  frame  is  used,  with  x  and  y  the  cross-shore  and 
alongshore  directions,  respectively.  Additional  pressure 
gages  in  8  m  depth  (x  -  800  m)  are  not  shown,  (bottom) 
Alongshore-averaged  seafloor  elevation  relative  to  mean  sea 
level  versus  x  (collected  during  August,  September,  and 
November). 

where  h  is  the  depth  and  g  is  the  gravity.  The  corresponding 
energy  £  and  cross-shore  energy  fluxes  T  of  shoreward  and 
seaward  propagating  waves  (superscripts  ±,  respectively)  at 
frequency/ and  location  x  are 

£*(/»  =  \  [< Copp(f,x )  +  (h/g)Co„(f,x). 

±  (2  \/h/gj  Copu(f,x)  J ,  (2) 

=  £±(f,x)y/ghi  (3) 

where  Copu  is  the  p  -  u  cospectrum  and  Copp  and  Co^  are 
p  and  u  autospectra,  respectively.  Unlike  the  individual 
shoreward  and  seaward  fluxes  (equation  (3)),  the  net  cross¬ 
shore  flux  (T+  -  T~  =  hCoPu )  does  not  require  the 
assumption  of  near  shore  normal  propagation  [Stoker,  1947]. 

[7]  The  pressure  and  current  meter  data,  sampled  at  2  Hz 
and  processed  in  3  hour  segments,  were  quadratically 
detrended  and  then  divided  into  448  s  demeaned  ensembles 
with  50%  overlap.  After  tapering  each  ensemble  with  a 
Hanning  window,  cross  spectra  and  spectra  with  about  48 
degrees  of  freedom  and  frequency  resolutions  of  0.002  Hz 
were  calculated.  Energy  flux  densities  (equation  (3))  were 
integrated  over  the  infragravity  frequency  band  (0.004- 
0.05  Hz)  to  estimate  bulk  infragravity  fluxes  F±  and  bulk 
reflection  coefficients  R2: 

0.05Hz 

F±(x)  =  J  F±(f,x)df,  R2(x)=F~(x)/F+(x)  .  (4) 

0.004Hz 


(i) 
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Figure  2.  (a)  and  (b)  Pressure  spectral  density  Copp\  (c)  and  (d)  normalized  phase  difference  0/tt 
between  pressure  at  jc  =  240  m  and  pressure  at  other  cross-shore  locations;  (e)  and  (f )  cross-shore  velocity 
spectral  density  Couu ;  (g)  and  (h)  normalized  phase  difference  0/tt  between  collocated  pressure  and  cross¬ 
shore  velocity;  and  (i)  and  (j)  infragravity  reflection  coefficient  R2  versus  cross-shore  location  (left)  7 
November  1997,  0700-1000  EST,  offshore  significant  sea-swell  wave  height  Hsig$m  =  2.0  m.  (right)  10 
September  1997,  0100-0400  EST,  Hsigy 8m  =  L0  m.  The  frequency  is  0.014  Hz.  The  circles  are 
observations,  and  the  solid  curves  are  solutions  of  the  linear  shallow-water  equations  for  a  normally 
incident  cross-shore  standing  wave  (R2  -  1)  on  the  measured  bathymetry.  The  predicted  Copp  and  Couu 
are  normalized  with  the  maximum  observed  value. 


Bulk  infragravity  energies,  fluxes,  and  reflection  coeffi¬ 
cients  were  obtained  from  each  collocated  p  -  u  sensor  pair. 
Estimates  from  the  same  alongshore  array  (Figure  1)  were 
averaged.  Errors  in  the  bulk  fluxes  and  R  owing  to  the 
assumption  of  normal  incidence  are  estimated  to  be  <20% 
(Appendix  A).  To  eliminate  observations  with  significant 
nongravity  motions  (e.g.,  shear  waves),  3  hour  runs  with  a 
mean  alongshore  velocity  >50  cm  s_1  at  any  current  meter 
were  excluded.  Frequency-alongshore  wave  number  spec¬ 
tra  (T.  J.  Noyes,  personal  communication,  2000)  confirm 
that  nongravity  wave  motions  contribute  generally  <20%  of 
the  infragravity  velocity  variance  in  the  680  runs  retained  in 
this  study.  Linear  finite  depth  theory  was  applied  to  pressure 
data  in  8  m  depth  to  estimate  Ess# m,  the  sea-swell  energy 
integrated  over  the  frequency  band  0.05-0.24  Hz.  Sig¬ 
nificant  wave  heights  in  8  m  depth,  tfSig,8m  =  4E]si%m,  in  the 
retained  3  hour  runs  ranged  from  0.2  to  2.7  m. 

[8]  Quadratic  difference-frequency  interactions  between 
sea-swell  components  with  slightly  different  frequencies 
(/  and  /  +  A/)  are  always  nonresonant  in  intermediate  and 
deep  (for  the  sea  and  swell  waves)  water  depths  and, 
theoretically,  result  in  a  bound  infragravity  wave  of  fre¬ 
quency  A f  that  does  not  satisfy  the  linear  dispersion 
relation.  In  this  nonresonant  case  the  doubly  integrated, 
normalized  bispectrum  (bn  in  equation  (9)  of  Herbers  et  al 
[1994]  and  p  here)  is  proportional  to  the  fraction  of  the 
total  infragravity  wave  energy  contained  in  bound  compo¬ 
nents.  In  contrast,  when  the  water  depth  is  shallow  for 
the  sea  and  swell,  near-resonant  quadratic  difference-fre¬ 
quency  interactions  can  result  in  gradual  (over  several 


wavelengths)  alteration  of  the  phase  and  energy  of  free 
infragravity  waves.  In  this  case  the  real  and  imaginary 
parts  of  the  bispectrum  are  related  to  the  rate  of  nonlinear 
phase  change  (wave  number  shift)  and  nonlinear  energy 
exchange,  respectively  [Herbers  and  Burton ,  1997;  Herb¬ 
ers  et  al,  2002].  Therefore,  in  the  shallow  depths  consid¬ 
ered  here,  p  (incorporating  both  real  and  imaginary  parts 
of  the  bispectrum)  is  interpreted  qualitatively  as  a  non- 
dimensional  measure  of  the  strength  of  near-resonant 
forcing  of  free  inffagravity  waves  by  quadratic  interactions 
with  sea  and  swell. 

[9]  Calculating  p  from  observed  pressure  or  velocity  time 
series  mixes  information  about  seaward  and  shoreward 
propagating  infragravity  waves.  Following  Elgar  and  Guza 
[1985],  nonlinear  coupling  also  was  estimated  from  time 
series  of  pressure  fluctuations  containing  only  shoreward  or 
only  seaward  propagating  infragravity  waves  (p+  and  p~, 
calculated  using  r\+  and  r\~  (equation  (1)),  respectively).  An 
error  in  accounting  for  the  effect  of  time  domain  windowing 
in  previously  reported  bispectral  integrals  [Herbers  et  al. 
1994,  1995a;  Ruessink,  1998]  is  corrected  below,  resulting 
in  p,  p+,  and  p“  values  that  are  a  factor  of  0.58  lower  than 
uncorrected  values. 

3.  Results 

[10]  Similar  to  previous  observations  [e.g.,  Suhayda , 
1974],  at  a  fixed  infragravity  frequency  the  observed 
cross-shore  variation  of  spectral  levels  of  pressure  (Figures 
2a  and  2b)  and  cross-shore  velocity  (Figures  2e  and  2f)  have 
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Figure  3.  (a)  and  (b)  Bulk  seaward  and  shoreward  sea-swell  energy  fluxes;  (c)  and  (d)  bulk  seaward  and 
shoreward  infragravity  energy  fluxes;  (e)  and  (f )  bulk  infragravity  reflection  coefficient  R2;  and  (g)  and 
(h)  measures  of  nonlinear  phase  coupling  p*  versus  cross-shore  location  x:(left)  7  November  1997, 
0700-1000  EST,  //Sig,8m  =  2.0  m  and  (right)  10  September  1997,  0100-0400  EST,  Hsi^ 8m  =  1.0  m.  The 
p+  biphases,  not  shown,  are  close  to  180°,  as  observed  by  Ruessink  [1998].  Circles  and  crosses 
correspond  to  shoreward  and  seaward  infragravity  wave  propagation,  respectively.  Note  the  different 
vertical  scales  for  fluxes  in  the  left  and  right  panels. 


minima  and  maxima  close  to  the  theoretical  locations  of 
standing  wave  nodes  and  antinodes. 

[u]  Note  that  the  observed  u  spectral  maxima  and  p 
spectral  minima  occur  at  the  same  location  (e.g.,  x  —  185  m), 
in  contrast  to  a  progressive  wave  field  where  the  spatial 
variations  of  u  and  p  spectral  levels  are  similar  and  do  not 
have  a  structure  with  nodes  and  antinodes.  The  observed 
phase  differences  between  spatially  separated  p  time  series 
(Figures  2c  and  2d)  and  between  collocated  p  and  u  time 
series  (Figures  2g  and  2h)  also  suggest  substantial  shoreline 
reflection.  With  a  progressive  wave  field  on  a  gently  sloping 
bottom  the  phase  difference  between  spatially  separated 
time  series  increases  approximately  linearly  with  increasing 
separation  along  the  propagation  direction,  and  the  phase 
between  collocated  p  and  u  is  zero.  At  the  representative 
infragravity  frequency  (0.014  Hz)  shown  in  Figure  2,  strong 
reflection  is  confirmed  by  the  R2(f)  estimates  (based  on 
collocated  pressure  and  cross-shore  velocity  observations 
equations  (2)  and  (3))  in  the  range  0.2-0.8  (Figures  2i  and 
2j).  As  discussed  by  A.  Sheremet  et  al.  (manuscript  in 
preparation,  2002),  similar  reflection  coefficient  estimates 
are  obtained  when  spatially  separated  sensors  are  included 
in  the  analysis.  Individual  frequency  bands  are  not  consid¬ 
ered  further  here. 

[12]  In  both  wide  (tfSig,8m  =  2.0  m)  and  narrow  (Hs^Sm  = 
1.0  m)  surf  zone  cases,  infragravity  wave  reflection  is 
significant  (Figures  3e  and  3f),  whereas  sea-swell  reflection 
is  weak  (Figures  3a  and  3b).  The  bulk  shoreward  infra¬ 
gravity  flux  F+  increases  shoreward  in  the  region  where  p+ 


is  large  (roughly  x  >  385  m  in  Figures  3c  and  3g  and  x  > 
175  m  in  Figures  3d  and  3h).  Within  the  surf  zone, 
decreases  to  negligible  values,  and  the  growth  of  F+  ceases 
(with  the  exception  of  one  observation  at  x  =  175  m  in 
Figure  3c).  Seaward  propagating  infragravity  waves  are  not 
coupled  to  groups  of  sea  and  swell  (pT  is  small  everywhere; 
Figures  3g  and  3h).  The  bulk  reflection  coefficient  R2 
(equation  (4))  decreases  from  values  as  large  as  1.5  far 
seaward  of  the  surf  zone  (Figure  3f)  to  between  about 
0.5  and  1.0  within  the  surf  zone  (Figures  3e  and  3f). 
At  the  location  nearest  the  shoreline,  R2  «  1.  These 
trends  are  consistent  with  the  hypotheses  [e.g.,  Longuet- 
Higgins  and  Stewart ,  1962]  that  F+  increases  in  the 
shoaling  region  owing  to  nonlinear  forcing  by  sea-swell 
groups,  that  sea-swell  breaking  reduces  the  nonlinear 
forcing,  and  that  shoreward  propagating  infragravity 
waves  are  reflected  strongly  at  the  shoreline. 

[13]  At  the  most  offshore  p-u  location  (x  =  500  m)  the 
seaward  F~  and  shoreward  F+  fluxes  both  increase  \yith 
increasing  Ess# m  (Figures  4a  and  4c).  For  all  Ess, 8m,  F  at 
jc  =  500  m  (5-6  m  depth)  and  x  =  210  m  (2-3  m  depth)  are 
approximately  equal  (their  ratio  is  about  1.2;  Figure  4b).  In 
contrast,  the  ratio  of  F+  at  x  -  2 1 0  m  to  F+  at  x  =  500  m  is  as 
large  as  3  with  moderate  Ess$ m  and  is  <1  with  the  largest 
Ess, sm  (Figure  4d). 

[14]  At  the  deepest  (5-6  m  depth)  location  the  depend¬ 
ence  of  the  bulk  reflection  coefficient  R2  on  ESSj 8m 
(Figure  5a)  is  similar  to  that  observed  in  13  m  depth 
[Herbers  et  al. ,  1995b,  Figure  5].  For  moderate  values  of 
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Figure  4.  (a)  Bulk  seaward  infragravity  energy  flux  F~  at  x  =  500  m;  (b)  bulk  ratio  of  F~  at  x  =  210  to 
K~  at  x  =  500  m;  (c)  bulk  shoreward  infragravity  energy  flux  F+  at  x  =  500  m;  and  (d)  bulk  ratio  of  F*  at 
x  =  210  to  F+  at  x  =  500  m  versus  sea-swell  energy  in  8  m  depth,  Ess#m. 


Ess, 8m>  R2  >  I?  whereas  for  the  largest  and  smallest  ESSi 8m, 
Br<  1 .  However,  the  maximum  R2  values  here  are  about  1 .5, 
lower  than  those  observed  (as  high  as  4.0)  in  13  m  depth. 
Closer  to  shore  (x  =  2 1 0  m),  R2  <  1  for  all  Ess# m  (Figure  5b). 
Although  values  are  scattered,  at  both  locations  the  mean  R2 


decreases  with  large  and  increasing  Ess$m.  At  both  locations 
the  nonlinear  coupling  \x+  increases  as  Ess# m  increases  from 
the  lowest  levels,  reaches  a  maximum  at  intermediate  Ess# m 
values,  and  then  decreases  with  further  increases  in  Ess$ m 
(Figures  5c  and  5d). 


x  =  500  m  x  =  210  m 


Figure  5.  (a)  and  (b)  Bulk  infragravity  reflection  coefficient  R2  and  (c)  and  (d)  bulk  nonlinear  phase¬ 
coupling  [i+  versus  sea-swell  energy  in  8  m  depth  ESSy$m  at  cross-shore  locations  (left)  x  =  500  and  (right) 
x  -  210  m. 
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Figure  6-  (a)  Average  values  of  bulk  shoreward  F+ 
(circles)  and  seaward  F~  (crosses)  infragravity  energy 
fluxes  normalized  by  Fx=5 oo?  the  seaward  infragravity  flux  at 
x  =  500  m,  and  (b)  bulk  nonlinear  phase  coupling  to  sea- 
swell  groups  of  the  total  (p,  triangles),  shoreward  propagat¬ 
ing  (p+,  circles),  and  seaward  propagating  (pT,  crosses) 
infragravity  wave  fields  versus  the  shoreward  sea-swell 
energy  flux  Fj  normalized  by  the  sea-swell  flux  in  8  m 
depth  Fss,Hm  •  The  vertical  bars  indicate  ±1  standard 
deviation  about  the  mean.  When  offshore  sea  and  swell 
energy  are  low,  the  dominant  infragravity  waves  often  are 
arrivals  from  remote  sources  [Herbers  et  al,  1995b].  To 
focus  on  locally  generated  infragravity  waves,  385  3  hour 
runs  with  Ess>% m  >  400  cm2  are  used. 


[15]  On  the  basis  of  similar  variations  of  observed  p,  with 
Ess ,  Ruessink  [1998]  suggested  that  the  intense  wave  break¬ 
ing  in  the  surf  zone  reduces  the  forcing  of  infragravity 
waves.  The  dependence  of  the  normalized  infragravity 
fluxes  and  nonlinear  coupling  on  the  reduction  of  sea-swell 
energy  flux  in  the  surf  zone  (F*  normalized  by  the  flux  in 
8  m  depth,  F+#m)  is  shown  in  Figure  6.  Values  of  the 
normalized  sea-swell  flux  F^s/Fss$m  close  to  1  correspond 
to  observations  seaward  of  the  surf  zone,  and  lower  values 
are  progressively  farther  within  the  surf  zone.  The  normal¬ 
ized  shoreward  infragravity  flux  F^iF*  =500m  increases 
as  F^/F*, 8m  decreases  from  1,  reaches  a  maximum  for 
F'ssIFss, 8m  ~  0.5,  and  then  decreases  in  the  inner  surf  zone 
(Figure  6a).  This  decrease  suggests  significant  damping  of 
infragravity  waves  in  the  surf  zone.  Similar  to  the  nonlinear 
coupling  of  the  total  infragravity  field  p  to  groups  of  sea  and 
swell  [Ruessink,  1998],  the  coupling  of  shoreward  propa¬ 
gating  infragravity  waves  p+  is  reduced  greatly  in  the  surf 
zone  (Figure  6b),  consistent  with  the  cessation  of  nonlinear 
forcing. 

[  16]  Although  the  dependence  of  infragravity  wave  prop¬ 
erties  on  the  normalized  sea-swell  flux  F+s  /Fss#m  demon¬ 
strates  the  important  role  of  sea-swell  breaking,  all 
observations  seaward  of  the  surf  zone  are  collapsed  into 
F&fFss£m  ~  1*  Consequently,  the  variation  of  infragravity 


fluxes  observed  in  the  shoaling  region  (i.e.,  x  >  175  m  in 
Figure  3d)  is  obscured.  To  include  this  spatial  variation, 
infragravity  wave  properties  are  shown  (Figure  7)  as  func¬ 
tions  of  the  normalized  cross-shore  position  x/x50 ,  where  x50 
is  the  location  where  F&  is  50%  of  F^,Sm  and  the  shoreward 
infragravity  flux  is  maximum  (Figure  7a).  F+/F'=50om 
increases  in  the  shoaling  region,  attains  a  maximum  at 
about  xix5o  ~  1*0,  and  then  decreases  (Figure  7a).  The 
coupling  p+  is  nearly  constant  through  the  shoaling  region, 
then  starting  at  the  outer  edge  of  the  surf  zone  (x/x50  «  2) 
decreases  sharply  to  negligible  values  in  the  inner  surf  zone 
(jc/jc50  <  1)  (Figure  7b).  Seaward  propagating  infragravity 
waves  are  not  phase  coupled  to  sea-swell  groups  (p_  always 
is  small;  Figure  6b),  and  the  energy  flux  of  the  unforced 
seaward  propagating  waves  F~  is  less  spatially  variable 
than  the  flux  F+  of  the  nonlinearly  forced  shoreward 
propagating  waves  (Figure  7a).  Dissipation  and  refractive 
trapping,  discussed  by  A.  Sheremet  et  al.  (manuscript  in 
preparation,  2002),  both  contribute  to  the  decrease  in  F~ 
with  increasing  offshore  distance  from  the  shoreline. 

[17]  The  dependence  of  R2  estimates  on  x/x50  is  shown 
in  Figure  8,  including  additional  estimates  farther  from 
shore  based  on  array  measurements  in  8  m  depth. 
Although  the  analysis  methods  are  different,  R2  estimates 
based  on  nearshore  collocated  p ,  u  observations  and  off¬ 
shore  array  measurements  (analyzed  without  the  assump¬ 
tion  of  normal  incidence  [Herbers  et  al,  1995b])  overlap 
smoothly  at  intermediate  jc/x50  values  (Figure  8).  Far 
offshore  of  the  surf  zone  (largest  x/x50),  F~  >  F+  and 
the  bulk  reflection  coefficient  R2  approaches  1.6  (Figure  8). 


Figure  7.  (a)  Average  values  of  bulk  shoreward  F+ 
(circles)  and  seaward  F~  (crosses)  infragravity  energy 
fluxes  normalized  by  F^=500,  the  seaward  infragravity  flux 
at  x  =  500  m,  and  (b)  bulk  nonlinear  phase  coupling  to  sea- 
swell  groups  of  shoreward  propagating  infragravity  waves 
(p+,  circles)  versus  x/jc50,  where  x50  is  the  cross-shore 
location  where  the  shoreward  sea-swell  flux  is  reduced  to 
50%  of  the  flux  in  8  m  depth.  The  dashed  vertical  line 
indicates  xlx50  =  1.  The  vertical  bars  indicate  ±1  standard 
deviation  about  the  mean.  To  focus  on  locally  generated 
infragravity  waves,  385  3  hour  runs  with  Ess#m  >  400  cm2 
are  used. 
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Figure  8.  Bulk  infragravity  reflection  coefficient  R2  versus 
normalized  cross-shore  location  jc/x50*  Circles  are  estimates 
based  on  collocated p  -  u  sensors  (equation  (1),  depths  <6  m), 
and  squares  are  based  on  pressure  array  observations  in  8  m 
depth.  The  385  3  hour  runs  with  Ess$ m  >  400  cm2  are  used. 

During  shoaling,  F+  increases,  but  F~  remains  relatively 
constant  (Figure  7a),  and  R 2  decreases  to  a  minimum  of 
about  0.5  near  x/x50  =  1.  In  the  inner  surf  zone  (x/x50  <  1) 
where  F+  decreases  toward  shore  (presumably  owing  to 
dissipation),  R2  increases  again  to  a  maximum  value  of 
about  0.8  at  the  shoreline  where  the  infragravity  waves  are 
strongly  reflected. 

f is]  Numerical  and  laboratory  studies  suggest  that  infra¬ 
gravity  energy  also  may  be  generated  by  fluctuations  of  surf 
zone  (Figure  7a)  width  and  setup  at  infragravity  periods 
[Symonds  et  al ,  1982;  Lippmann  et  al. ,  1997].  Although 
such  generation  is  not  precluded  by  the  results  presented 
here,  much  of  the  observed  growth  in  F+  occurs  well 
seaward  of  the  surf  zone  and  thus  offshore  of  the  region 
where  generation  by  a  time-varying  breakpoint  theoretically 
occurs. 

4.  Directionality  of  Infragravity  Waves 

[19]  The  method  used  here  to  estimate  the  shoreward  and 
seaward  propagating  infragravity  wave  components  is  based 
on  the  assumption  of  near  shore  normal  propagation. 
Although  the  energy  flux  and  R2  estimates  are  not  degraded 
significantly  by  the  infragravity  wave  obliquity  in  these  data 
(Appendix  A),  obliquity  could  be  associated  with  refractively 
trapped  edge  waves.  As  in  previous  studies  in  similar  water 
depths  [e.g.,  Huntley  et  al.,  1981;  Oltman-Shay  and  Howd , 
1993],  frequency  and  frequency- wave  number  spectra  show 
deviations  from  normal  incidence.  In  cases  with  large  mean 
incident  angles  the  spectral  levels  of  v  often  are  within  a 
factor  of  2-3  of  the  cross-shore  velocity  u  levels  (Figure  9a). 
At  infragravity  frequencies,  frequency -alongshore  wave 
number  spectra  of  p ,  u,  and  v  (p  is  shown  in  Figure  9b)  have 
variance  maxima  at  wave  numbers  within  the  edge  wave 
range.  On  the  basis  of  frequency-alongshore  wave  number 
spectra  and  the  qualitative  standing  wave  character  of  the 
cross-shore  variation  of  infragravity  phases  and  variance 
(Figures  2b,  2d,  2f,  and  2h),  observations  similar  to  these 
have  been  interpreted  as  consistent  with  edge  waves,  the 


normal  mode  solutions  of  gravity  waves  undergoing  multiple 
reflections  between  the  shoreline  and  an  offshore  turning 
point.  However,  the  cross-shore  analysis  for  this  run  shows  a 
substantial  progressive  shoreward  energy  flux  (Figure  3f) 
that  suggests  that  a  significant  fraction  of  the  energy  is  lost 
before  waves  reach  the  shoreline.  This  motivates  the  consid¬ 
eration  of  strongly  damped  edge  waves  by  A.  Sheremet  et  al. 
(manuscript  in  preparation,  2002). 

5.  Conclusions 

[20]  Observations  of  pressure  and  velocity  collected  in 
water  depths  between  1  and  6  m  are  used  to  estimate  the 
seaward  and  shoreward  energy  fluxes  of  infragravity  waves. 
Although  there  is  considerable  scatter,  on  average  the 
shoreward  infragravity  flux  increases  during  shoaling  when 
the  nonlinear  coupling  between  shoreward  propagating 
infragravity  waves  and  sea-swell  groups  is  strongest.  Within 
the  surf  zone,  coupling  is  reduced  strongly,  and  the  shore¬ 
ward  flux  decreases,  consistent  with  the  cessation  of  forcing 
and  increased  dissipation.  Reflection  is  strong  at  the  shore¬ 
line.  The  coupling  of  seaward  propagating  infragravity 
waves  to  sea-swell  groups  always  is  weak.  The  results 
support  the  existing  hypothesis  that  shoreward  propagating 
infragravity  waves  are  amplified  in  the  shoaling  zone  by 
nonlinear  interactions  with  sea-swell  groups  and  are 
strongly  reflected  near  the  shoreline.  Damping  of  infragrav¬ 
ity  waves  in  the  surf  zone  also  appears  to  be  important.  The 
seaward  infragravity  flux  decreases  with  increasing  distance 
from  the  shoreline,  consistent  with  both  dissipation  and 
refractive  trapping,  but  is  less  spatially  variable  than  the 
shoreward  flux.  As  a  result  of  these  spatially  varying  energy 
fluxes,  the  bulk  infragravity  reflection  coefficient  exceeds  1 
far  seaward  of  the  surf  zone,  reaches  a  minimum  of  about 
0.4  in  the  outer  surf  zone,  and  increases  to  about  0.8  at  the 
shoreline  where  infragravity  waves  are  reflected  strongly. 

Appendix  A:  Error  Estimates  for  R2  and  F* 

[21]  Errors  in  infragravity  energy  fluxes  and  reflection 
coefficients  estimated  using  collocated  observations  of 
pressure  p  and  cross-shore  velocity  u  and  the  assumption 
of  normal  incidence  (equations  (2)  and  (3))  are  examined 
using  linear  WKB  theory  for  shallow-water  waves.  The 
cross-shore  and  the  alongshore  coordinates  are  x  and  y ,  and 
the  x  origin  is  at  the  shoreline.  Assuming  that  the  depth 
varies  only  in  the  cross-shore  direction,  h  ~  h(x ),  the 
pressure  time  series  can  be  expressed  as  a  superposition 
of  Fourier  modes  with  different  frequencies  /  and  along¬ 
shore  wave  numbers  ky: 

00  00 

p(t,x,y)  =  J  df  J  P'{f,x,ky)^y^ 

OO  OO 

or,  equivalently,  incidence  angles  9 (/  x): 

OO  it 

P(t,x,y)  =  J  dfe-W  J  P(f,x,Q)dk^dQ,  (Al) 

— OO  — Tt 

where  P  is  the  cross-shore  structure  of  a  Fourier  mode  and 
k  —  k(f  x)  is  the  local  wave  number  modulus.  Normal 
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ic  (1/m)  p2  (m2s) 

Figure  9.  (a)  Frequency  spectra  of  infragravity  cross-shore  u  (circles)  and  alongshore  v  (triangles) 
velocities,  (b)  normalized  frequency-cyclic  alongshore  wave  number  (f  -  k)  spectra  of  pressure  p 
estimated  with  the  Maximum  Entropy  Method  [Wu,  1997],  and  (c)  frequency  spectra  of  pressure  at  x  = 
310  m,  on  10  September,  0100-0400  EST.  The  dashed  curves  in  Figure  9b  are  edge  wave  dispersion 
curves  for  modes  0-5  (approximately  the  boundary  between  the  leaky  and  trapped  wave  domains).  The 
/-  k,  spectra  at  each/are  normalized  by  the  frequency  spectra,  and  darker  shades  indicate  higher  spectral 
levels  within  each  frequency  band.  Additional  results  for  this  run  are  shown  in  Figures  2(right)  and 
3(right). 


incidence  corresponds  to  0  =  0.  A  relation  similar  to 
equation  (Al)  can  be  written  for  U,  the  cross-shore  structure 
of  the  cross-shore  velocity.  P  and  U  are  assumed  to  have  the 
form 


/>(/>,  9)  = 


+  j, 


U(f,x, 0)  =  ] cos  0, 


(A2) 


where  A±  =  A\f  x,  0)  are  the  real  amplitudes  of  shoreward 
(+)  and  seaward  (— )  propagating  waves,  e±  is  an  arbitrary 
phase,  and  <j>  is  the  phase  with  respect  to  the  shoreline 


(A3) 


where  0  =  0(/  x)  is  the  incidence  angle  at  the  seaward  edge 
of  the  integration  domain  and  kx  is  the  cross-shore  wave 
number,  respectively.  Shoreward  propagating  waves  are 
assumed  to  be  partially  reflected  at  the  shoreline,  with  the 
WKB  phases  €+  —  €_  =  —it/ 2  determined  by  matching  to 
exact  standing  wave  solutions  [Herbers  et  al.,  1995a]. 
Using  equation  (A2),  the  autospectra  and  cross-spectra  for 
collocated  pressure  and  cross-shore  velocity  time  series  are 

n/2 

XBP(/,x)=i  J  [(/T)2  +  (A-  f+2AJ  A-  sin(26)p0, 

-*/2 


*/2 

X„„(/,x)  =  l-  J  [(^+)2  +  (A-f  -  2A+A-  sin(24>)] cos20d0, 
-*/2 
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Figure  Al.  Contours  of  (a)  Rq,  (b)  the  magnitude  of  the  maximum  value  attained  by  |(3|  over  a  plane 
beach  with  a  0.02  slope,  over  the  depth  range  1  -6  m,  and  (c)  Fq/F,+  as  functions  of  the  mean  direction 
0m  and  directional  spread  a  (angular  spread  function  5(0)  is  defined  in  equation  (A  14)).  Each  dot 
represents  the  alongshore-averaged  bulk  infragravity  mean  direction  and  directional  spread  for  a  single 
3  hour  run  at  a  given  cross-shore  location,  computed  with  a  method  that  utilizes  both  components  of 
horizontal  velocity  [Herbers  et  al.,  1999],  (d)  The  relative  errors  of  R2  and  F±  for  R%  =  0.005,  |(3|  = 
0.0015,  and  FZ/Ft+  =  1.1  (e.g.,  the  relative  error  for  R2  is  defined  as  (R2-R,2)/R2)  versus  the  true 
reflection  coefficient  R,2. 


„(/»  J  [{A+f  -  (A-)2  +  2^+yf_cos(24>)] coseae 


[22]  Substituting  equation  (A4)  into  equations  (2)~(3) 
yields  the  estimated  energy  fluxes  at  frequency  f: 


ir/2 

:F±=^r  j  {[(^+)2+(^“)2](1+cos20) 

-V2 

±2  [(^4+)2  -  (4-)2]  cos  0  +  2^+/f-(24>)sin2e|(/0.  (A5) 

The  amplitudes  are  related  to  the  energy  £ +  and  angular 
distribution  S  of  the  shoreward  propagating  modes  by 

{A+f  =  2£+{f,x)S{f,S\ 

(A~)2  =  2R]£+  {f  ,x)S(f ,%),  (A6) 

A+A~  =2Rl£+(f,x)S(f,Q), 

where  JZ^2  5(/,  0)  d()  —  1 .  Specular  reflection  and  a  true 
reflection  coefficient  Rt  that  is  independent  of  frequency  and 
direction  are  assumed  for  simplicity.  With  £+  and  S 


assumed  independent  of  frequency  over  the  infragravity 
band  the  frequency-integrated  fluxes  are 


F±  =  \Z-E 


Uv2 

[(1  -f  i?2)(l  +  cos20)  ±  2(1  —  i?2)cos0j 

it/2  } 

E  J  df  J  5(0)5,  sin(2(J>)  sin2  0rf6  > 


Df  —  tt/2 


where  E+  is  the  total  shoreward  propagating  energy  in  the 
infragravity  band  Df  of  width  A/. 

[23]  To  assess  the  errors  in  equation  (A7),  it  is  convenient 
to  normalize  by  the  true  shoreward  energy  flux 


F+  =  yfghE+  J  5(0)  cos  QdO  , 


yielding 


n2  r>2\  F«  / 


j+  =  (1  +  Wt  +  Rffi)  j£(=  1  if  F+  is  exact), 
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F~ 

Ff 

— 


,F+ 


(R20  + 13 R,  +  R2)  (=  R]  if  F  is  exact) 


1  +3^r‘  +k2rI 

1+3 R,+RJRl 


(=  1  if  R2  is  exact) , 


(A9) 

(A10) 


where  R2,  and  Ff  are  the  values  of  the  estimated  reflection 
coefficient  and  shoreward  flux  for  R,2  =  0  and 


^°+b)-  (A,,) 


with 


tv/2  ti/2 

i=  J  S(0)(1  +  cos29)d0,  b=  J  25(0) cos 

— tt/2  -h/2 

V2 

•  =  2A/  j  df  J  5(0)  sin  (24>)  sin2  0d0. 


0d0 


(A12) 


D,  —tx/2 


Assuming  that  the  beach  slope  is  constant,  the  integral  (A3) 
becomes 


HM)  = 


I'kJx 

Vsh(x)  l- 


COS0  4-  -r-r 

sm0 


(A13) 


[24]  The  errors  in  the  estimates  (A8)-(A10)  can  be 
calculated  for  a  given  Rt2  if  values  of  Rq,  (3  and  F^IF^ 
are  known.  These  depend  on  £(0)  and  are  obtained  by 
substituting  in  equation  (A  12)  an  analytic  expression  for  S 
[Longuet-Higgins  et  al ,  1963], 

r  ”i  2(2/az— 1) 

S(0)  =  A  [cos(0  -  0W ) j  ,  (A  1 4) 

where  A  is  a  normalization  constant,  0m  is  the  mean 
propagation  direction,  and  a  is  the  directional  width.  The 
dependence  of  Rq,  (3,  and  Fq/F+  in  equations  (A8)-(A10) 
on  0^  and  a  is  shown  in  Figures  Ala-Alc.  The  infragravity 
band  was  taken  as  0.01-0.03  Hz  (a  wider  bandwidth 
reduces  the  errors),  and  the  beach  slope  was  taken  as  0.02 
(approximately  the  average  Duck  slope  over  the  span  of  the 
sensor  array).  Figure  Alb  shows  the  magnitude  of  the 
maximum  of  |(3|  over  the  depth  range  1-6  m.  For  typical 
observed  0W  and  a,  values  of  |(3|,  Ro,  and  F$/Ft+  are  about 
0.005,  0.0015,  and  1.1,  respectively.  With  these  values  the 
relative  errors  for  R2  and  F*  do  not  exceed  20%  for  0.05  < 
Rt 2  <  4  (Figure  Aid).  Note  that  because  the  fluxes  F±  are 
overestimated  by  approximately  the  same  amount,  the  flux 
ratio  R2  is  more  accurate  than  the  fluxes 
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